Sensing distant nuclear spins with a single electron spin 
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We experimentally demonstrate the use of a single electronic spin to measure the quantum dynam- 
ics of distant individual nuclear spins from within a surrounding spin bath. Our technique exploits 
coherent control of the electron spin, allowing us to isolate and monitor nuclear spins weakly coupled 
to the electron spin. Specifically, we detect the evolution of distant individual ^^C nuclear spins 
coupled to single nitrogen vacancy centers in a diamond lattice with hyperfine couplings down to 
a factor of 8 below the electronic spin bare dephasing rate. Potential applications to nanoscale 
magnetic resonance imaging and quantum information processing are discussed. 



Detection and control of the magnetic signals from 
nuclear spins is an important problem in science and 
technology. Single nuclear spin detection remains an 
outstanding goal in magnetic resonance imaging (MRI), 
and could have far-reaching implications from physics 
to medicine [1-3]. Likewise, nuclear spins stand out in 
quantum information science for their exceptional iso- 
lation from their environment, making them attractive 
qubit candidates. However, as this isolation results from 
their weak magnetic moments, individual nuclear spins 
are extremely difficult to detect and control. Utilizing the 
electronic spin associated with a single nitrogen- vacancy 
(NV) center in diamond to both detect and control sur- 
rounding nuclei is a promising approach to this challenge, 
and significant efforts are currently directed towards this 
goal. Early work demonstrated that NV centers can be 
used to sense strongly coupled proximal nuclear spins [4- 
6]. Subsequently, this approach has been used to cre- 
ate a few-qubit quantum register [7, 8], perform simple 
quantum algorithms [4, 9], implement single shot read- 
out of both nuclear and electronic spins [10, 11], and 
demonstrate a multi-second quantum memory at room 
temperature using a single coupled ^^C nuclear spin in 
isotopically purified diamond [12]. 

In this Letter we show that the NV electronic spin 
can be used to isolate and probe the quantum dynamics 
of distant, weakly coupled individual Carbon- 13 nuclear 
spins. Our approach relies on dynamical decoupling se- 
quences, which enhance the sensitivity to individual nu- 
clear spins while suppressing NV electronic spin deco- 
herence [13-19]. Crucially, this allows us to observe the 
coherent evolution of nuclear spins whose coupling to the 
NV is weaker than the limit imposed by the NV spin's in- 
homogenous dephasing rate, l/T^ . We demonstrate this 
capability by identifying and measuring the coupling to 
single nuclear spins from amongst a bath of naturally 
abundant ^^C nuclear spins. The ability to isolate spins 
from within a bath is important for applications where 
surrounding nuclear spins cannot be avoided; moreover, 
it enables the use of these nuclear spins as a resource. In 
particular, our technique could be extended to harness 
weakly coupled nuclear spins in quantum registers [4, 8- 
12], to investigate the spatial extent of the NV electronic 



wave function [20-23], or to detect individual spins out- 
side of the diamond lattice for single spin MRI applica- 
tions [1, 2, 13, 24]. 

The central idea of this work is depicted in Fig. Ic, and 
can be understood in terms of the coherent evolution of a 
single ^^C nuclear spin interacting with the NV electronic 
spin sensor. Through their interaction, the magnitude 




FIG. 1: (a) Schematic setup. We measure shallow- implanted 
NV centers weakly coupled to ^^C nuclear spins. A dc mag- 
netic field is applied using a movable permanent magnet, (b) 
Pulsed spin manipulation and readout, (c) Central concept. 
Conditional evolution of the ^^C on the Bloch sphere dur- 
ing NV spin manipulation, for ^^C spin initially in state W). 
Upper (lower) panel shows ^^C evolution with the NV spin 
initially in state |0) (|1)). Gray boxes show evolution during 
spin echo (single vr-pulse). Additional 7r-pulses push condi- 
tional nuclear spin evolution further apart, increasing its en- 
tanglement with the NV and thus the signal. For clarity a 
sequence with an odd number of vr-pulses is shown; the result 
is qualitatively the same for the even numbered sequences 
used in this work, but is more difficult to visualize. 
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and orientation of tfie local magnetic field experienced 
by the ^^C spin depends on the NV spin state. As a re- 
sult, when the NV spin is prepared in a superposition of 
its energy eigenstates, the ^^C spin undergoes conditional 
Larmor precession around two different axes. During this 
precession the electron and nuclear spin become entan- 
gled and disentangled, resulting in collapses and revivals 
in NV spin coherence [5]. However, if the NV- nuclear 
coupling is weak then the two precession axes are similar 
and the resulting entanglement is small. The key idea 
of this work is to increase the degree of entanglement, 
and thus the measurable signal, by applying periodic tt- 
pulses, flipping the NV spin with a frequency matched to 
the precession of the ^^C spin. As shown schematically 
in Fig. Ic, this constructively enhances the conditional 
evolution of the ^^C state further apart, increasing its 
entanglement with the NV spin. 

To describe the experimental sequence shown in 
Fig. lb, we consider the mutual interaction between a 
single NV electronic spin-1 and a single ^^C nuclear spin- 
i, which is governed by the Hamiltonian {h = 1) 



H = ASl - iieB.S, + (m„B + S,A) ■ I. 



(1) 



Here, S (I) and /ig (Mu) are the electronic (nuclear) 
spin and magnetic moment respectively, B is the exter- 
nal magnetic field, and A is the hyperfine interaction. 
Due to the large zero field splitting A/27r c:^ 2.87 GHz, 
we have made the secular approximation conserving Sz- 
The nuclear spin evolves conditionally on the NV spin 
state according to i^nuc[^s] = ^^^^^ • n^^, where a 
is the vector of Pauli spin matrices. Thus the nu- 
clear spin precesses about an effective magnetic field 
axis n^^ = (/inB + rrisA) /uorns with Larmor frequency 
^rus = l/inB + m^AI when the NV spin is in state Im^). 
We initialize the NV spin in a superposition of |0) and 
|1) and apply a periodic series of tt pulses spaced by 
evolution time t, during which |0) and |1) accumulate a 
relative phase. Converting the phase into a population 
difference, the normalized fluorescence signal at the end 
of the measurement of total length A^t is p = (5 + l)/2, 
where 



5 = 1 - (no X ni)2 sin^ (^) sin^ (^) 



sin^ (Ar^/2) 
cos2 (0/2) 



(2) 



and cos (j) = cos ^ cos ^ — no • ni sin ^ sin This 
is the extension for arbitrary (even) pulses of the well- 
known result for = 1 corresponding to Hahn spin echo 
[5]. The factor of A^ multiplying the small angle (j) is 
responsible for the enhanced signal that enables detection 
of weakly coupled nuclear spins. We use Eq. (2) with 
A^ = 8 to fit all data below. 

Experiment. — We study NVs implanted in a diamond 
sample with naturally abundant spin-^ ^^C (1.1% abun- 
dance). The NV electronic orbital ground state is a spin 
triplet which can be initialized in the \ms = 0) state us- 
ing green illumination, and read out via spin-dependent 



fluorescence [5]. The \ms = il) degeneracy is lifted by 
a dc fleld applied using a permanent magnet. We coher- 
ently manipulate the NV spin using resonant microwaves 
delivered by a coplanar waveguide. The control pulse se- 
quence used in this work is known as XY4-8 [14], and 
consists of A^=8 7r-pulses as shown in Fig. lb, phase al- 
ternated to mitigate microwave pulse errors. 

Detection of a single ^^C spin with a hyperflne coupling 
of |A| /27r ~ 125 kHz is shown in Fig. 2a. The dip in 
the NV fluorescence signal at evolution time t ^ 9 fis 
results from the coherent interaction with a single ^^C 
spin and represents the signal. Importantly, this signal 
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FIG. 2: Coherent interaction of an NV-Center with a single 
^^C detected using XY4-8. (a) The evolution time between 
TT-pulses t is swept with a fixed magnetic field of 15 G, data 
(black) and fit (red). The blue line is the calculated response 
without the individual ^^C. The inset shows the calculated 
response to Hahn echo, XY4-4, and XY4-8, with the collapse 
resulting from the background spin bath kept fixed at the 
calculated response to XY4-8 to highlight the impact of the 
single ^^C. (b) The fluorescence is plotted as a function of t 
and magnetic field strength B, with magnetic field oriented 
along the NV z axis (a = 0, as defined in the inset). The 
dashed line highlights the data shown in (a), (c) Fluorescence 
as a function of t and magnetic field orientation a, with B 
= 19 G. The black line displays the extracted dip position in 
both (b) and (c). Here the total coupling strength is | A| /27r ^ 
125 kHz. 
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would go undetected using spin echo, as shown in the 
inset of Fig. 2a. In addition to the dip arising from the 
single ^^C, an overall decay is apparent and arises from 
the interaction with the bath of distant ^^C spins. The 
impact of this spin bath can be treated as an effective 
fluctuating ac magnetic field at the characteristic Larmor 
precession frequency [5] . Consequently, the NV evolution 
will also periodically exhibit collapses and revivals in its 
coherence due to the bath ^^C spins, and the visible decay 
in Fig. 2a is the onset of the first collapse (for the exact 
form of the collapse used in the fits, see [25]). 

To confirm that the observed feature arises from a sin- 
gle ^^C, we plot the measured signal as a function of 
applied magnetic field strength (Fig. 2b) and orientation 
(Fig. 2c). We measure the NV response in both possible 
sets of magnetic sublevels (|0) ^ |1) and |0) ^ | — 1), 
see [25]) and the entire data set is simultaneously fitted 
to Eq. (2). The only free parameters are the hyperfine 
vector A, defined relative to the NV axis and magnetic 
rotation axis (see [25]), and two free parameters describ- 
ing the strength and spectral width of the surrounding 
spin bath [25]. The calculated signal using the extracted 
fit parameters is shown in the right panels, and is in good 
agreement with the data. The fit gives a total coupling 
strength of |A| /27r ~ 125 kHz, which is less than the 
bare dephasing rate of this NV, l/T^ = 400±16 kHz. 
T2* is measured using a Ramsey sequence with no control 
TT-pulses (measurement shown in [25]), and arises from 
the bath spin configuration varying from measurement to 
measurement. This demonstrates that we can simultane- 
ously decouple the NV from the ^^C spin bath and still 
obtain a measurable signal from a target nuclear spin. 

Figure 3 shows the measurement of the weakest cou- 
pled individual ^^C observed in this work. This single 
spin produces no discernible signal for the experimental 
settings in Fig. 2 (measurement shown in [25]), because 
the weaker hyperfine coupling results in a dip at longer 
evolution times which occurs after the first collapse 
induced by the bath and is therefore suppressed. In or- 
der to measure this weakly coupled spin at very long 
evolution times we increase the magnetic field, moving 
a bath-induced revival into the range of interest. We 
thereby observe a characteristic dip in the bath revival 
resulting from an individual nuclear spin. Again, this 
data is reproduced using our model calculations; the ex- 
tracted coupling strength is |A| /27r 47kHz, a factor of 
- 4.6 below the measured l/T^ = 217 ± 9 kHz [25]. 

The coherent nature of the hyperfine interaction also 
allows us to observe the simultaneous impact of multi- 
ple ^^C spins coupled to a single NV. As the maximum 
mutual interaction strength of two ^^C spins on adjacent 
lattice sites is 2 kHz [26] , it is negligible compared to the 
hyperfine interaction and to leading order we can treat 
the -^^C spins as independent. Consequently, the signal 
given by Eq. (2) becomes a product of the individual con- 
tributions from each coupled spin, S = HjSj [5]. Figure 4 
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FIG. 3: Detection of a distant ^^C with a XY4-8 sequence. 
As in Fig. 2 the NV fluorescence is shown as a function of (a) 
magnetic field strength B with a = 0.26 rad, and (b) magnetic 
field orientation a, with B = 80 G. The dip from the single 
^^C is visible within a spin-bath revival centered at 2-2t{/ujq. 
The extracted coupling strength is |A| /27r 47kHz. 



shows the impact of three ^^C spins on the evolution of 
an NV center. The double peak visible in the data re- 
sults from two overlapping dips caused by two ^^C spins 
with similar coupling strengths. The peak in the middle 
reflects the product of two negative contributions and di- 
rectly reflects the coherent nature of the interaction. This 
is emphasized in the inset of Fig. 4c, where the calculated 
effect of all three individual ^^C spins on S is displayed 
separately on top of the bare signal resulting from the 
bath. The weakest coupled of the three ^^C spins has an 
extracted coupling strength of |A| /27r ^ 64kHz, a factor 
of - 8 below the measured l/T^ = 560 ± 60 kHz [25]. 

Discussion. — We now turn to a discussion of our re- 
sults and the limits of our technique. In contrast to the 
nearby ^^C spins detected in previous works [5, 22, 23], 
the magnitudes of the hyperfine couplings observed here 
are consistent with a purely dipolar interaction. Assum- 
ing a point-dipole interaction, the measured couplings 
translate into NV-^^C distances of ~ 0.4-0.8 nm [25], the 
range expected for a probabilistic distribution of natu- 
rally abundant ^^C. However, ab initio calculations sug- 
gest that at this length scale the spatial extent of the 
electronic wavefunction can lead to significant contribu- 
tions from the contact interaction [20-23] . Because our 
observed coupling strengths are an order of magnitude 
weaker than the contact interaction strengths predicted 
by these calculations [20, 21], we conclude that the dipo- 
lar interaction provides a minimum distance, and that the 
observed ^^C spins may be farther away than estimated. 
These considerations imply that our technique provides a 
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and (b) NV luminescence is shown vs. evolution time t and 
magnetic field orientation a (data and fit), the applied field 
strength is 19 G. The fit assumes three ^^C spins and yields 
unique individual coupling strengths, while fits assuming two 
or fewer ^"^C spins could not reproduce the data. The NV re- 
sponse along the dashed line is shown in (c) black, red: data, 
fit. The traces in the inset display the calculated response of 
the NV interacting separately with each of the three coupled 
^^C spins, and with the spin bath alone (black.) 



tool to measure the spatial distribution of the electronic 
wave function beyond the limits of current theory, al- 
though a larger data set in the spirit of recent surveys 
of strongly coupled ^^C spins [22, 23] would be required 
to provide the statistics necessary to identify individual 
distant ^^C lattice sites. 

The sensitivity limits of our technique are determined 
by the NV spin decoherence. The hyperfine couplings 
|A| /27r extracted in this work exceed the decoherence 
rate I/T2 obtained in a Hahn spin echo sequence, and in 
principle the ^^C spins could therefore be detected using 
only a single 7r-pulse [25]. However, this would require 
precise optimization of the magnetic field strength and 
orientation, which is not practical in a typical experi- 
ment. One advantage of extended pulse sequences is to 
relax these conditions and greatly simplify the identifica- 
tion of nuclear spins whose hyperfine coupling with the 
NV spin are initially unknown. As a result, using the pre- 
sented technique we were able to detect and characterize 
the closest ^^C for every NV investigated. 

The use of dynamical decoupling sequences with mul- 
tiple TT-pulses has the added benefit of decoupling the 
NV spin from the surrounding environment, increasing 
its coherence time from T2 to Tf^ = T2N'^^^ and thereby 



potentially improving the sensitivity beyond the limit set 
by T2. Increasing the number of 7r-pulses for the same 
total evolution time requires higher magnetic fields to 
ensure that the sequence remains synchronized with the 
evolution of the target nuclear spin. We note that in 
the limit | /in B| |A| a weakly coupled single ^^C has 
the same Larmor precession frequency as the surround- 
ing nuclear spin bath and cannot be isolated. Therefore 
I A| /27r > I/T2 is the limit for the detection of individual 
^^C nuclear spins in a sample of natural isotopic abun- 
dance. However, for applications involving the detection 
of spins with a different gyromagnetic ratio than the ^^C 
bath we find that the sensitivity limits of our technique 
improve with number of 7r-pulses as N^^^ [25]. As a re- 
sult, our approach could potentially be used to detect 
individual spins inside or outside the diamond lattice up 
to the ultimate limit |A| /27r > l/2Ti, where 1/Ti is the 
NV spin relaxation rate. At room temperature l/2Ti 
can be as low as ~100 Hz [27, 28], corresponding to the 
dipolar coupling between an NV and a proton spin at a 
distance of ~7-9 nm, which is of great interest for single 
spin MRI applications [1-3]. 

In conclusion, we used coherently controlled single NV 
electronic spins to detect distant nuclear spins. Extended 
pulse sequences allow for detection of individual weakly 
coupled ^^C spins that would otherwise be unresolvable. 
We also showed that the simultaneous detection of sev- 
eral distant ^^C spins is possible, even within an envi- 
ronment consisting of a large number of spins with the 
same gyromagnetic ratio. Our technique allows for sensi- 
tive, coherent measurements of the nuclear spin environ- 
ment of a single electronic spin; moreover, extensions of 
this approach could be used to exploit the nuclear spin 
environment as a resource to extend the size of control- 
lable multi-spin quantum registers. Potential novel appli- 
cations range from information storage to environment- 
assisted sensing and single spin MRI. 

Authors' Note. — Following submission of this work two 
complementary studies appeared [29, 30] in which the 
technique demonstrated in this work was used to resolve 
individual ^^C nuclei weakly coupled to a single NV in 
the high magnetic field regime [29], and in isotopically 
purified diamond with depleted ^^C nuclei [30]. 
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SYSTEM HAMILTONIAN AND SIGNAL DERIVATION 



Because of the large zero-field splitting of the NV-spin S and moderate external magnetic fields, its interaction with 
the nearby nuclear spin I is treated within the secular approximation (S = Szz). The effect of the external magnetic 
field on the electron is undone by the echo sequence and therefore neglected {h= 1). The reduced Hamiltonian is 

iy = 5^A-I + /inI-B (1) 

As the Hamiltonian is diagonal with respect to the NV spin, the nuclear spin evolves conditionally on the NV spin 
state in an effective magnetic field B^^ = dJ^^n^j {cums = 

n-B -|- TTT/g A| , Urns 

Hnnclrris] = ^msl"^^ ' (2) 

In the experiment, the NV spin states = ±1 are detuned, we can therefore treat the NV spin as a two-level system; 
here we restrict ourselves to = and vfis = 1. In the applied echo sequence the NV spin state is initially prepared 
in the direction (in the rotating frame); the nuclear state is in a completely mixed state, the initial state of the 
system is therefore given by the density matrix. 

Po = \{^s^(yx)®\^i (3) 

As the final 7r/2-pulse rotates about the same axis as the initial one, we measure the expectation value (cr^) (referred 
to as S in the main text). The fiuorescence rate (experimental signal p) is simply derived from this as 

F=^(l + (0) (4) 
The expectation value is rewritten as (5R designates the real part) 

(a,) = Tt[pa,] = 5RTr[p{t)a+] = 25RTr„,c (0|p|l) (5) 
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This is structurally similar to an interference amplitude of the two NV spin states |0) , subject to some evolution 
described by p{t). Throughout the sequence, the NV spin is flipped = ^ = 1; in total N times, N is even 
(odd) in case of XY4 (PDD). Figure SI schematically shows the microwave manipulation of the sequences. Depending 
on the sequence the rotation axes are either x (PDD and CP) or alternating x and y (XY4). The time evolution is 
generally given by 



U = . . .exp[-^Hnuc[ms]r2] exp[-i7r/2cr3,] exp[-^i^nuc[^s]n] 

Inserting (3), (6) in (5) and utilizing the commutator relations 

ex.-p[i7T/ 2a x,y]H nuclei = i^nuc[l] exp[i7r/2cra;^^] and vice versa 
ex.-p[i7T/2(jy](Jx = — cTaj exp[i7r/ 2 cr^y] 



(6) 



(7) 



we obtain 

(<T^) = l^{0\{ts + (-1)V^)|1) Tw[exp[iiJnuc[0]T„] exp[iiJ„uc[l]r„-i][. . .] exp[-iiJ„uc[0]r„_i] exp[-iiJnuc[l]T„] 



(8) 



where j designates the total number of applied 7r-pulses around the y-axis. This is the only difference between a 
Carl-Purcell (CP-A^) and XY4-A^ sequence [1]; in accordance with most publications we ignore that sign and set it to 
+1. The expression will now be evaluated in the case of PDD-A^ and (XY4)-A^ sequences. We identify repeating parts 
of the sequence and describe their total evolution by a single effective rotation around n with angle (/), both depending 
on the evolution time t. From the last equation it is clear that a XY4-A^ sequence can be split into smaller repeating 
parts (referred to as XY below) (see flg. Sib). In the experiment the XY4 sequence is utilized, as it decreases the 
impact of pulse imperfections compared to XY or PDD. 

The repetition of these parts can then conveniently be carried by multiplying the total angle by the number of 
repetitions k. In case of XY4 (PDD) it applies k = N/2 {k = {N ^ l)/2) 



(^rep)^ = (l COS + j 



1 COS k(j) + ian sin k(j) 



(9) 



In equation (8) the product of two such rotations is regarded, arising from the |0) and |1) state of the NV spin. With 
the assumption that both rotations employ the same effective angle 0, yet different rotation axes n, m we obtain (in 
the last step the Chebyshev polynomials are utilized). 



UfJJrh =l(cos^ k(j) + nm sin^ k(j)) + it7((n + m) sin k<p cos k(j) — n x m sin^ k(j)) 

=IL(T|(cos(/)) + nmU^_-^{cos (I))) + i^^((n + m)T/e(cos (/))/7/e_i(cos (/)) — n x mU^_-^{cos (j))) 

Inserting this into eq. (8) we obtain (utilizing Pell's equation T|(cos(/)) + cos^ (j)U^ _^{cos (p) = 1) 

(ax) =T^(cos0) — nmUl_i{cos <p) = 1 — (1 + nm — cos^ (j))Ul_i{cos c 
nm 



=T2(cos< 



(l-Tfc^(cos< 



1 — cos 

The second equation is useful to determine the boundaries of the signal for all possible number of repetitions k. 



(10) 

(11) 
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FIG. SI. Schematic Sequences, (a) XY4-A^ sequence with N being the total number of vr-pulses. As the rotation axes x and y 
have the same effect on the time evolution operator, it can be deduced that the periodicity is shorter as indicated in (b). (c) 
shows the PDD sequence. 
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XY4 Sequence 



In the case of XY4 the iterative part of the sequence can be expressed as (see fig. SI); we use the abbreviation 

0m, = ^mtl2 



and 



U\ \t I 2]Uq [t]Ui [t 1 2] =1 (cos 00 cos 0i — no • ni sin 0o sin 0i ) 

+it7 • (—no sin 0o cos 0i — fii sin 0i cos 0o + 2ni x (lii x no) sin^ 0i/2 sin 0o) 
=1 cos 6 + icrn sin = 1 cos + icrn 



Z//q [t/2]Z//|[t]Z//Q [t/2] =1 (cos 00 cos 01 — no • ni sin0o sin0i) 

+icr • (fio sin 0o cos 0i + lii sin 0i cos 0o — 2no x (lio x lii ) sin^ 0o /2 sin 0i ) 
=1 cos + icrm sin 



(13) 



(14) 



Inserting n, m and cos0 into eq. (11), we obtain the exphcit result with cos0 given by the preceding equations. 

(cr^) = 1 - 4(1 - (noni)2)sinVo/2sinVi/2(l -cos0)/7^_i(cos0) (15) 

We now regard the solution written in the form (12). As < < 1] < 1, the fraction in the equation represents 

the lower limit of the signal independent of number of repetitions /c, see fig.S2. In particular, the absolute minimum 
is obtained when the axes of the effective rotations n, m are parallel. As they have the same length, they are identical 
n = m. Inserting this into eqs. (14), (15) yields the two conditions. 



sin 00 cos 01 + sin^ 0i /2 sin 0o — noni sin^ 0o/2 sin 0i = 
sin 01 cos 00 + sin^ 0o /2 sin 0i — fioni sin^ 0i /2 sin 0o = 



(16) 
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FIG. S2. Signal strengths for varying number of evolutions (green, red, blue: XY 2,4,6; respectively); the black curve is the 
lower boundary to all the possible XY sequences. 



PDD Sequence 

The evolution can be written as 

V(o[t]V(i[t] =1 (cos 00 cos 01 — no • ni sin0o sin0i) (17) 
-\-ia • (no sin 0o cos 0i + lii sin 0i cos 0o — sin 0o sin 0ino x lii) 
=1 cos + icrn sin = 1 cos + icrn 

m is obtained by reversing the angles 0o —00, 0i — 01? inserting this into eq. (11) yields the result for PDD. 

(a^) = 1 - 2(1 - (noni)2) sin' 0o sin' 0i/7^'_i(cos 0) (18) 

Which is easily seen to reduce to the well known result for Hahn spin-echo as Uk-i^o = 1- 
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BACKGROUND DUE TO THE SPIN BATH 



We treat the impact of the bath of ^^C as a classical magnetic field, whose frequency distribution S{uj) is Gaussian 
around the Larmor precession due to the external field, giving rise to decaying collapses and revivals. Please note that 
this does not include the dominant part of the decoherence, which is a result of low- frequency noise due to fiip-fiop 
processes within in the bath, see e.g. [2], [3]. The coherence envelope ((cr+)) in the system may be written as 



((^+)) exp[-x[t]] exp 



du F{ut) 



27r cj2 

Here F{z) defines the sequence-dependent filter function; it takes the form 
F{z) = ao + ^ an cos(nz/2) 



(19) 



(20) 



Assuming the width a of the Gaussian distribution to be narrow a <C cjo we can simplify it as 

S{uj) = ~X^ujI exp[-(cj - cjl) V(2^')] - ^ A'cj' exp[-(a; - uji^f / {2ct^)] 
x{t) can be easily integrated piecewise, yielding 

x{t) = ^ao + ^ an exp(|n^t^cr^) cos(ntcjo)^ 
The high frequency components in the filter function die out faster, as shown in fig. S3. 



(21) 



(22) 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 
Evolution time l/w 

FIG. S3. Simulated bath decoherence, XY-2,4,8 (black, blue, red, respectively) a = 0.1c<;/(27r), A = 0.25. Faster oscillations 
can be seen to be damped more quickly. 

In addition, the low frequency components in magnetic noise lead to a dominant overall decay in coherence. Increas- 
ing the number N of 7r-pulses in a sequence reduces the susceptibility to low frequency noise and thereby enhances 
the effective coherence time Tf^, see. [1],[4]. 



((^+))lowfreq=exp 



exp 



(t|«) 



(23) 



We include this term when we regard a revival of the bath induced coherence (fig. 3 in the main part) and use it for 
scaling arguments when analyzing the limits of our method. 



APPROXIMATION OF THE HYPERFINE INTERACTION 



For the following discussion, we regard the situation in which the coordinate system is chosen such that the z- 
axis coincides with the symmetry axis of the NV center, the NV spin's center of mass defines the origin, and the 
^^C lies within the x-z plane, (a special case of the more general case shown in fig. S4). Motivated by its small 
magnitude, we treat the hyperfine coupling as rotationally symmetric about reducing the hyperfine interaction A 
to A = {Ax^O^ Az}. Note that the contact interaction is still treated exactly within this simplification. As the NV 
wavefunction is not perfectly rotationally symmetric, it can potentially produce a magnetic field at the ^^C location 
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that has components in the y direction, which we neglect here. From the consistency of our fits we infer that these 
terms are smah. 

Previous works have observed an enhanced ^^C gyromagnetic ratio for transverse fields, arising from coupling to 
the NV electronic spin (see for example [5].) This enhancement of the gyromagnetic ratio is a result of corrections to 
the secular approximation that we do not consider as they are negligible when the coupling A is weak enough such 
that /ie//^n <^ ^/ is the case for all ^^C spins considered in this work. 




FIG. S4. Definition of the coordinate system. The z-axis is defined by the symmetry axis of the NV center, while the x-axis is 
defined such that the magnetic field vector B always lies within the x-z plane as the magnet is rotated. The angle of magnetic 
field rotation a is defined relative to the z-axis, but the arc is shown here without the first 7r/2 segment from the z to x axes 
for clarity. The position of the ^^C relative to the NV is then given by the spherical coordinates r, 0, and 0, while A is the 
hyperfine coupling vector resulting from this spatial configuration. 



SIGNAL STRENGTH VS. NUMBER OF APPLIED PULSES 



When the hyperfine coupling strength |A| exceeds the bare Hahn echo decoherence rate I/T2, a strong signal can 
be obtained without enhancing the coherence to Tl^ using extended sequences. We therefore numerically investigate 
the maximal impact Pmax of an individual ^^C in the regime I/T2 > |A| using extended sequences. 

At first we will consider the impact arising form the coherent evolution described by eq. (15) which is suppressed 
by low frequency noise, see eq. (23). We therefore regard an expression of the form 

P=Hl-(<^.)](W)) (24) 

We calculate the maximum signal Pmax for a given coupling strength |A|, the orientation and strength of the external 
magnetic field as well as the evolution time t are optimized numerically. Figure S5 shows the result. The signal 




Number of repetitions k 



FIG. S5. Numerically calculated signal strength vs. number of pulses and coherence time; the dashed lines highlight the scaling 

as k^/' 

strength can be seen to scale with the number of repetitions k as k^^^; this can be understood as follows: ideally, the 
two states of the ^^C superposition get pushed apart linearly with the number of repetitions k and the evolution time 
t. For small differences in these states, this leads to a quadratic increase in signal strength. One therefore obtains 



p (xt^k^ exp 



dp 

m 



= ^ tock^^^ and Pniax oc k'^^^ 



(25) 
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If the investigated spin is a ^^C, one has to include the cohapses and revivals in coherence coming from the Larmor 
precession (at frequency cjl of the ^^C bath spins. For cjl ^ |A| the signal p tends to peak att~7r/co'L(l/2 + n), n: 
integer; however the signal in this range is exponentially suppressed due to the collapses from the bath. 

RAMSEY MEASUREMENTS 

We infer the inhomogeneous dephasing rate I/T2 of the NV spin from Ramsey measurements, as shown in fig. S6. 
The result can be accurately reproduced taking only the coupling to the nuclear spin of the nitrogen into account 
(coupling strength: 2.2 MHz). We measure up to an evolution time of 40 ps to ensure that there are no revivals in 
coherence. 




Evolution time [[js] Evolution time [\is\ Evolution time [[js] 

FIG. S6. Ramsey measurements of the NV centers investigated in the main text. Data shown in (a), (b) and (c), corresponds 
to the NV data shown in fig. 2, 3 4 in the main text. Extracted T2 are 2.5 ± 0.1, 4.6 ± 0.2, 1.8 ± 0.2 ]is, respectively. 



CONSISTENCY OF RAMSEY AND XY4 MEASUREMENTS 

As a consistency check of our method, we show in fig. S7 measurements performed on an NV center with a strongly 
coupled i^C, where A > l/T^. We measure its impact on the NV spin evolution in a XY4-8 measurement (a,b) as 
well as in a Ramsey sequence (c,d). The extracted hyperfine couplings are 422 kHz (Ramsey) and 438 kHz (XY4-8); 
these values are larger the inhomogeneous dephasing rate 286 kHz. 
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FIG. S7. XY4-8 and Ramsey measurements, (a) and (b) display XY4-8 measurement and fit vs. magnetic field orientation and 
strength, respectively, (c) shows a Ramsey measurement performed on the same NV, it is fitted incorporating an additional 
energy splitting due to the coupling to the ^^C. (d) is the fourier transform of the data shown in (c). In addition to the 2.2 
MHz splitting from the nitrogen, a clear splitting due to the hyperfine coupling is visible; lines mark the position of the peaks 
as determined from (c). 



ADDITIONAL XY4-8 MEASUREMENTS - COMPLEMENTING FIGS. 2,4 



In fig. S8 we show the remaining data from the data sets used to generate the fits shown in figs. 2 and 4 of the main 
text. In general, we measure the NV response against magnetic field strength and orientation for both the upper and 
lower NV spin transitions. 
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FIG. S8. Measurements and simulation completing the data sets shown in the main text, showing the NV response vs magnetic 
field strength (a,c,d) and orientation (b). (a and b) Signal when addressing the lower transition of the NV shown in fig. 2 (main 
text). (c,d) The data set corresponding to the NV shown in fig. 4 of the main text; (c): upper transition, (d): lower transition. 



8 



ADDITIONAL XY4 MEASUREMENTS - COMPLEMENTING FIG. 3 

In fig. S9 we show additional data for the NV with the weakest coupled ^^C observed, (fig. 3 of the main text). The 
experimental settings are similar to those of figs. 2 and 4 in the main text, but XY4-16 is used in place of XY4-8 to 
further enhance any impact of the weakly coupled ^^C. No clear impact from the ^^C is observed, which is consistent 
with our calculated response based on the fit displayed in the main text. 
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FIG. S9. (a) Measurement of NV subject to a XY4-16 sequence, (b) Calculated response (solid lines) of the NV obtained by 
fitting (22) to the data shown in (a); dashed lines include the additional effect of the isolated nearby ^^C based on the coupling 
parameters extracted from the fit to the data shown in fig. 3 in the main text. 



PARAMETER TABLE OF INVESTIGATED NV CENTERS 

In table I we present the extracted parameters from the NV centers we observed that could be treated within the 
framework derived above. In particular, NVs with strongly coupled ^^Cs (|A| > IMHz) were not investigated in 
detail. The hyperfine coupling strengths A are determined from the fits to the XY-4 measurements; the extracted 
distance r and angle 6 (defined in fig. S4,) are based on the assumption of point-like dipole interactions, see main text 
for details and discussion. 

NV index (arb.) T2 [ps] number of ^^Cs [kHz] [kHz] |A| [kHz] distance r [A] angle [rad] appears in fi g. 



1 2.5 1 100 -75 126 5.6 18 2,S6,S8 

2 4.6 1 -40 -22 46 8.0 -22 3,S6 

3 3.5 1 -222 -359 422 3.6 -10 S7 

4 2.0 2 180 -51 187 5.1 27 

63 -61 88 6.3 15 

5 1.8 3 98 97 138 6.3 60 4,S6,S8 

115 -13 116 6.1 32 4,S6,S8 

37 -52 64 6.9 12 4,S6,S8 



TABLE I. Extracted parameters from XY4-8 and Ramsey measurements. 
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